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ABSTRACT 

Flutter analysis of a two degree of freedom airfoil in compressible flow is performed 
using a state-space representation of the unsteady aerodynamic behavior. Indicial response 
functions are used to represent the normal force and moment response of the airfoil. The 
structural equations of motion of the airfoil with bending and torsional degrees of freedom 
are coupled to the unsteady airloads and the aeroelastic system so modelled is solved 
as an eigenvalue problem to determine the stability. The aeroelastic equations are also 
directly integrated with respect to time and the time-domain results compared with the 
results from the eigenanalysis. A good agreement is obtained. The derivatives of the 
flutter speed obtained from the eigenanalysis are calculated with respect to the mass and 
stiffness parameters by both analytical and finite-difference methods for various transonic 
Mach numbers. The experience gained from the two degree of freedom model is applied 
to study the sensitivity of the flutter response of a wing with respect to various shape 
parameters. The parameters being considered are: (i) aspect ratio, (ii) surface area of the 
wing, (m) taper ratio and (iv) sweep. The wing deflections are represented by Chebyshev 
polynomials. The compressible aerodynamic state-space model used for the airfoil section 
is extended to represent the unsteady aerodynamic forces on a generally laminated tapered 
skewed wing. The aeroelastic equations are solved as an eigenvalue problem to determine 
the flutter speed of the wing. The derivatives of the flutter speed with respect to the shape 
parameters are calculated by both analytical and finite difference methods. 
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1. INTRODUCTION 


Flutter is a dynamic aeroelastic instability that involves the coupling of inertial, elastic 
and aerodynamic forces. Studies on transonic flutter are important in the design of air- 
crafts which operate in the transonic regime. Moreover, as compressibility effects are very 
important in transonic flow, a compressible aerodynamic theory must be used to model 
the unsteady aerodynamic behavior. 

The various methods that axe used for flutter analysis of a two dimensional airfoil 
differ in the prediction of aerodynamic loads. The lift and moment predictions on an 
airfoil undergoing harmonic motion have been obtained by Theodorsen [1]. Several CFD 
methods which are used to determine the transonic flowfield around two dimensional airfoils 
are listed by Ballhaus and Bridgeman [2]. Ballhaus and Goorjian [3] performed time- 
marching transonic flutter predictions using the transonic aerodynamic code LTRAN2. 
Yang at al. [4] performed flutter analysis of the NACA 64A006 airfoil with pitching and 
plunging degrees of freedom using aerodynamic coefficients obtained from the transonic 
codes UTRANS2 and LTRAN2. Guruswamy and Yang [5] used the LTRAN2 code for 
studies on aeroelastic time response analysis of airfoils. The aerodynamic force and moment 
response of an airfoil can be represented by indicial functions. Indicial response for an 
incompressible flow was obtained theoretically by Wagner [6]. Jones [7] used a two-pole 
exponential approximation to the Wagner function. Venkatesan and JYiedmann [8] have 
given a three-pole indicial response function that can express the Theodorsen’s function 
over the entire reduced frequency range. 

Flutter analysis of wings are also done using different representations of the aerody- 
namic loads. Kapania, Bergen and Barthelemy [9] have used the Yates modified strip 
analysis [10] for representation of the unsteady aerodynamic loads in obtaining the flutter 
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response of a laminated wing, A recent review of various studies in unsteady transonic 
flow calculations was given by Tjatra [11], Advanced codes such as XTRAN3S [12] and 
CAP-TSD [13], which use the transonic small disturbance equation are currently being 
used for aeroelastic analysis. 

Leishman and Nguyen [14] have represented the aerodynamic indicia! response func- 
tions for compressible flow by upto three-pole approximations, the response consisting of 
two parts, one due to non-circulatory loading and the other due to circulatory loading. 
This has advantages over the CFD-based methods in the sense that the CFD methods axe 
in genera! computationally very expensive. 

In recent years, considerable efforts are being made to integrate the aerodynamic, 
structural and control aspects of the design of an aircraft. Livne [15] notes that for in- 
tegrated multidisciplinary wing synthesis, where design for aeroservoelastic stability is an 
objective, it is required to represent the aeroelastic equations of motion in Linear Time 
Invariant (LTI) state-space form. The unsteady aerodynamic loads on the wing can be 
represented in a state-space form, thereby adding only a small number of states to the 
mathematical model of the aeroservoelastic system. 

Sensitivity analysis is becoming an important design tool in engineering design ap- 
plications. It was first recognized as a useful tool for assessing the effects of changing 
parameters in mathematical models of control systems. The gradient based mathematical 
progra mm ing method used in optimal control and structural optimization furthered the 
development of sensitivity derivatives, because sensitivity derivatives are used in search 
directions to find opt imum solution [16]. Sensitivity analysis has also become a versatile 
design tool, rather than just an instrument of optimization programs [17]. Sobieski [18] 
discusses in detail about the System Design Derivatives which help in understanding the 
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effect a particular design variable would have on the desired performance of the system, if 
it were perturbed by a small percentage from its original value. 

The sensitivity derivatives of a system can be found using either analytical or finite 
difference methods. Analytical sensitivity analysis has found increased interest in engi- 
neering design as it eliminates uncertainity in the choice of step size needed in the finite 
difference method. The step size if too large leads to truncation errors and if too small 
leads to ill-conditioning. 

Adelman and Haftka [17] have shown that structural sensitivity analysis has been 
available for over two decades. Structural sensitivity analysis has been sufficient in the 
past because sizing variables such as plate thickness and cross-sectional areas affect the 
mass and stiffness properties of the airframe, but, not its basic geometry. Therefore, 
aerodynamic sensitivity analysis capability has been limited in development un til recently. 
For example, Rudisill and Bhatia [19] developed expressions for the analytical derivatives of 
the eigenvalues, reduced frequency and flutter speed with respect to structural parameters 
for use in minimizing the total mass. 

Pedersen and Seyraman [20], examined the change in flutter load as a function of 
change in stiffness, mass, boundary conditions or load distribution. They showed how 
sensitivity analysis can be performed without any new eigenvalue analysis. The solution 
to the main and an adjoint problem provide all the necessary information for evaluating 
sensitivities. Their paper ma i nl y focused on column and beam critical load distributions. 

Hawk and Bristow [21] developed aerodynamic sensitivity analysis capabilities in sub- 
critical compressible flow. They first analyzed a baseline configuration, and then calculated 
a matrix containing partial derivatives of the potential at each control point with respect 
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to each known geometric parameter by applying a first order expansion to the baseline 
configuration. The matrix of partial derivatives is used in each iteration cycle to analyze 
the perturbed geometry. However, this analysis only handles chordwise perturbation dis- 
tributions, such as changes in camber, thickness and twist. A new approach, which is 
still under development, has been presented by Yates [22] that considers general geometric 
variations, including planform, and subsonic, sonic and supersonic unsteady, nonplanar 
lifting-surface theory. 

Recently, Livne et al [23] applied an equivalent plate structural modeling, which in- 
cludes transverse shear, to an HSCT wing. Simple polynomials were used for Ritz functions 
and depth and thickness distributions. The derivatives of the stiffness and mass matrices 
were obtained analytically with respect to the shape variables of the wing. Livne [24] 
observed that as higher order polynomials are used for better modeling of the structure, 
the more sensitive is the finite difference derivative to the step-size used and in some cases, 
it is impossible to obtain any valuable information by finite differences. 

Barthelemy and Bergen [25] explored the analytical shape sensitivity derivatives of 
the wing’s aeroelastic characteristics, such as section lift, angle of attack, rolling moment, 
induced drag and divergence dynamic pressure, for subsonic subcritical flow, with respect to 
geometric parameters. Results showed the characteristics nonlinearity to be s mall enough 
to be well approximated by sensitivity based linear approximations. These approximations 
are valid within a range that is useful to designers in the initial design phase. 

Kapania [26] has obtained sensitivity derivatives of the flutter speed of a two dimen- 
sional airfoil in incompressible flow with respect to the mass and stiffness parameters. 
Kapania, Bergen and Barthelemy [9] have obtained the shape sensitivity derivatives of the 
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flutter response of a laminated wing in incompressible flow. In this work, Yates’ mod- 
ified strip analysis [10] was used for the aerodynamic model in conjunction with Giles’ 
equivalent plate analysis [27,28] for the structural model. 

Sensitivity derivatives are of great importance in integrated multidisciplinary design 
optimization of aircrafts. Karpel [29] used a gradient-based constrained optimization on a 
composite active-flexible wing to achieve aircraft performance requirements and sufficient 
flutter and control stability margins with a minimum weight penalty and without violating 
the design constraints. The sensitivity derivatives of the flutter dynamic pressure, control 
stability margins and control effectiveness with respect to structural and control design 
variables were obtained analytically. 

Hajela et al [30] applied Sobieski’s Global Sensitivity Equations (GSE) in an aircraft 
synthesis problem where the constraints involved the coupled disciplines of structures, aero- 
dynamics and flight mechanics. The coupled system was represented by smaller subsystems 
and the total behavior sensitivities were determined by applying the GSE method. 

Barthelemy et al [31] discuss a multidisciplinary design optimization method applied 
to a supersonic transport wing. Aerodynamic and structural disciplines are integrated for 
a minimum weight design under static aeroelastic constraints. He points out that as the 
number of dependent variables in each discipline becomes large, the calculation of the finite 
difference derivatives contributes substantially to the total optimization cost. 

In this paper, the flutter speed of a two degree of freedom airfoil with plunging and 
pitching degrees of freedom in transonic flow is determined. The aerodynamic force and 
moment response are represented by indicial response functions as given by Leishman and 
Crouse [32]. The resulting aerodynamic state equations are coupled with the structural 
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equations and the stability of the aeroelastic model so obtained is determined using an eige- 
nanalysis. Flutter calculations axe also performed in the time-domain using the Wilson-0 
method and the results compared with the eigenanalysis results. The flutter speed of the 
airfoil is estimated and the sensitivity derivatives of the flutter speed with respect to various 
parameters, namely, mass ratio, static unbalance, radius of gyration, bending frequency 
and torsional frequency axe calculated by both analytical and finite difference methods. 
The aerodynamic state-space model [14] is then modified to represent the unsteady aero- 
dynamic forces on a wing. The wing structure is modelled as a wing box and Chebyshev 
polynormals axe chosen for the displacement function. The aeroelastic equations for the 
wing axe solved as an eigenvalue problem to determine the stability. The derivatives of the 
flutter speed are calculated with respect to the shape parameters, namely (i) aspect ratio, 
(ii) area, (iii) taper ratio, and (iv) sweep, by analytical and finite difference methods. To 
the best of our knowledge, this is a first study on the sensitivity analysis of the flutter 
response in transonic flow. 
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2. AERODYNAMIC MODEL 


In recent yeaxs, considerable efforts are being made to integrate the aerodynamic, 
structural and control aspects of the design of an aircraft. Since the control and the 
structural dynamic behaviors can easily be expressed in terms of the state-space form (i.e., 
in terms of a set of first order ordinary differential equations in time), it is desirable that 
the unsteady aerodynamic airloads be also expressed in the same form. In recent yeaxs, 
considerable efforts have been made in that direction. 

The state-space approach has the advantage that any system of differential equations 
can be represented by a set of first order ordinary differential equations of the form 

x = Ax 4- Bu 

with the output equations given by 

y = Cx + Du ( 1 ) 

where x are the aerodynamic state variables, u axe the system inputs and y axe the system 
outputs. K the unsteady aerodynamic behavior can be represented by state equations, 
then they can be easily coupled to the structural equations of motion and the resulting 
system can be examined for aeroelastic stability. 

The aerodynamic force and moment response of an airfoil can be represented by indicia! 
functions. For example, Jones [7] used a two-pole exponential approximation to the Wagner 
function [6] given by 

<f>(S) = 1 - 0.165exp(— 0.04555) - 0.335exp(-0.35) (2) 

where S = 2Vt/c, V is the freestream velocity, t is the time and c is the chord. The 
state equations describing the unsteady aerodynamic response can then be obtained by the 
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application of Laplace transforms to these indicial functions. The resulting state equations 


are 
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4- 0.5a3 /^(t) 


( 4 ) 




where C n is the normal force coefficient and ct is the angle of attack. 


In this paper, the state-space representation given by Leishman and Nguyen [14] has 
been used to represent the compressible unsteady aerodynamics. The indicial normal force 
and quarter chord pitching moment responses to a step change in angle of attack a and a 
step change in pitch rate q can be written as [14] 


Cn(S) = ^i(S,M) + C N .(M)^(S,M) 


a 


“ = - M) + C Na (M)<^£(S, Af)(0.25 - x ac (M )) 


( 5 ) 


where are exponential functions of S and M. Here, M is the 

Mach number, q = occjV is the pitch rate, C n is the normal force coefficient, Cm is the 
pitching moment coefficient about the quarter chord and C Na is the normal force curve 
slope. The superscripts C and I refer to circulatory and non-circulatory components of 
the indicial response functions. Note that (where f3 = \/l - M 2 is the compressibility 
factor) in [14] has been replaced by C^ a (M), so that experimental values of C^ a obtained 
as functions of Mach number cam be used. 
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The aerodynamic state equations have been shown by Leishman and Nguyen [14] to 


be given by 


where 


x = Ax + B 


A — diag[a\\ 022 <*33 a 44 a ss a 66 a n a 8 s] 


1 110 110 0 

0.5 0.5 0 1 0 0 1 1 


The output equations are given by 


= Cx + D 


where 

Q _ Cll CJ2 C\z C14 0 0 0 0 

C21 C22 0 0 C25 C26 C27 C 28 

D f 4/M 1/M 

—1/M —7/12 M 

The nonzero terms of the dj j’s and Cjj’s are given in the Appendix. 
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3. FLUTTER ANALYSIS OF THE AIRFOIL 


AEROELASTIC MODEL 

The aerodynamic equations in state-space form can be coupled to the structural equa- 
tions of motion of an airfoil section with bending and torsional degrees of freedom. The 
equations of motion for the airfoil section shown in Fig.l can be written as 

mh + S&0 + g h h + rrujj 2 h h — Q h 

( 8 ) 

S$h + 1 $8 + god ■+- IocjqO = Qo 

where m = 7r/xp(c/2) 2 is the mass per unit length, p is the mass ratio, p is the air density, 
1 $ m(c/2) r| is the polar moment of inertia about the quarter chord per unit length, r 9 

is the radius of gyration about elastic axis, S$ = m(c/2)x^ is the static mass moment, xq 
is the nondimensional distance in semichords from elastic axis to center of mass, h is the 
plunge displacement (positive downward), 8 is the pitch angle, u> h and u>o are the bending 
and torsion frequency respectively, g h and g e are the structural damping coefficients in 
plunging and pitching respectively and Qk and Qo are generalized aerodyn ami c forces in 
plunging and pitching respectively. 

By defining the states 


zi = h, z 2 = 8 , z 3 = k, z 4 = 8 , 


the above equations can be written as 


where 


’I 0 ' 


' 0 1 


f° ) 

0 M 

Z = 

-lc - g 
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m 

So 
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( 10 ) 
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The resulting set of first-order differential equations in terms of the z and x states are 
given by 

i o oi r z ) r o i oifz'i 

0 M 0 < 1= Di-k Dj - g C' { \ (13) 

0 0 ij [x J [ B' Bj Aj (xj 

which is a 12x12 system of linear equations. The stability of the system could be determined 
at different free-stream speeds by an eigenanalysis of the above system of equations. The 
flutter speed is that particular value of the free-stream speed at which the real part of the 
eigenvalue approaches zero. 
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The aeroelastic equations (13) could also be integrated with respect to time using a 
time-integration scheme. The Wilson-0 method [33] was used for this purpose. A set of 
first order ODEs can be represented as 

[*]{»'} = [SIM (14) 

In the Wilson-0 method, it is assumed that the variation of acceleration from time t to 
t + 0A<, where 0 > 1.0, is linear. At time ( t + 0Ai), then 

0At 

{ut+e At } = {u t } + — Wt+ $At + u't} (15) 

Then equation (14) becomes 

[[*] - [51^] {« W = + *-£<} (16) 

Using the starting values of {u*} and {uj} at time is computed from equation 

(16). The vector is then calculated from equation (15). The new values of 

{^t+0A*} and are then used in equation (16) to update the vector. The 

step-by-step integration of the equations is done in this mann er with respect to time by 
repeating the above process. The amplitudes of plunge and pitch displacements Eire then 
monitored as time progresses. 
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4. FLUTTER ANALYSIS OF THE WING 


STRUCTURAL MODEL 

The structural formulation is based on a Ritz solution technique using the energy 
functionals for a laminated plate which includes the bending and stret chin g of the reference 
surface. The planform geometry can be represented by any generally tapered skewed 
configuration. The original rectangular (x, y) coordinate system and the transformed (77, £) 
coordinate system of the wing are shown in Fig.2. The x - y plane is the mid-plane of 
the wing and the z axis is normal to the wing. For an unswept wing the fiber angle is 
measured counterclockwise from the positive y axis. As the wing is swept, the fiber angle 
is also rotated correspondingly. 

In the Rayleigh-Ritz formulation, Chebyshev polynomials Ti are used to represent the 
displacements at any point on the wing[34]. The Chebyshev polynomials axe given by 

W) = 1 

= 0 (17) 

Ti(lp) — 2 ^i-2 — 1 < xj) < 1 

The displacements are expressed in terms of the Chebyshev polynomials as shown 

I J 

''(f.O'EEwira 0 

i=o j = 0 

KL 

(is) 

k —0 1=0 

M N 

W(v,0 = E P ™T m (v)T n (0 - 1 < < 1 

m= 0 n=0 

It has been shown by Singhvi and Kapania [34] that for free vibrations of the l amin ated 
composite wing (t.e., in the absence of aerodynamic forces) the equations of motion can 
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be derived using classical plate theory in the form 


[AflW + (*]{«} = 0 


( 19 ) 


where [K] and [M] axe the stiffness and mass matrices. The eigenvector {5} is defined as 


q — (Roo, R 01 , ■••Rij', Sqo, Soi, ... Ski‘,Poo * ^01 j •••P mn) 


( 20 ) 


Linear and rotational springs of large magnitude are placed at the wing root to satisfy the 
clamped boundary conditions. The stiffness matrix for the plate alone (i. e., excluding the 
springs) is 


m = £ Jjb) t [t] t 


A B 
B D 


[T][B]\J\d V d( 


( 21 ) 


where [ 5 ] is the matrix whose elements consist of the partial derivatives of the Chebyshev 
polynomials with respect to the natural coordinates 77 and f and is defined by 



( 22 ) 


where 



V v V( W v Wrfr, U>££ W 


The [T] in equation(21) is the transformation matrix that relates the strain and curvature 
vector in the (x, y) coordinate system to the strain and curvature vector in the (77, f) 
coordinate system and J is the Jacobian of the transformation. The strain transformation 
is given by 



( 23 ) 


where 



V>y "1“ XV 2ti)jy) 
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The details of the [T] and [5] matrices and J are given in [35]. A typical element of the 
mass matrix [M] is given by 

M >i = T i ( v )T l (0T4 n )T p (()IJ\J V d( (24) 

The coefficients R t] and Ski in {q} corresponding to the inplane displacements in 
equation(19) are condensed out using static condensation to the form 

{M]{P mn } + [K){P mn ) = 0 (25) 

where [. M ] is the mass matrix and [. K ] is the stiffness matrix of order (m + l)x(n + 1) with 

generalized coefficients {P mn }. In the present work, a value of 5 is chosen for both m and 
n. 

AEROELASTIC MODF.T. 

The aerodynamic state space model which was used for the aeroelastic analysis of a 
typical section is extended to represent the unsteady aerodynamic forces acting on a wing 
in transonic flow. The lift and moment forces on a typical section acting at the quarter 
chord are given by equation(lO) as 

L = ipK’cCV 

M = ipKVCW 

When extending this compressible aerodynamic theory to a finite span wing, the lift 
forces are assumed to be distributed along the quarter chord line (reference line) and the 
moments act about the reference fine. Since the lift and moment forces are non-conservative 
forces, using the principle of virtual work, we get 

= J —LShdy + J MSddy (26) 
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where l is the length of the quarter chord line, Sh and 86 are virtual displacements and y 
is the coordinate along the reference line. 

The displacement at any location y is given by 


h(y) = w(rj y £) (27) 

where 17 and £ are the natural coordinates corresponding to the (x,y) coordinates of the 
point at distance y from the origin. 

The rotation about the reference lin e is given by 

6(y) = w, x cosA - w, y sinA (28) 

For facilitating in numerical integration using Gaussian quadrature, the limits of integra- 
tion along the reference line are transformed in the range of —1 to 1 by y = /(I -f- 
where — 1 < if) < 1 


Substituting the expressions for the lift and moment on the wing and the wing deflec- 
tion, we have 

f LShdy 

Jo 

=ir 

2 j-i 

/ 1 m n r 


LShdip 

m n r /•! 


SR } 


and 


f MSddy 

Jo 

= L£ 

- L A^t W>’ ]{ 


MSOdip 

ll . . _ 2 L JL r r\ 


(29) 


x } +■ 




8P t , 
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where m and n represent the order of the Chebyshev polynomial used in the displacement 
function. The row vectors [C\ p ] and [C 2 P ] are the elements of the [C] matrix (see section on 
typical section) where p = 1, 2, . . . , 8. The row vectors [T>i p '] and [D 2p ,] axe the elements 
of the matrix given by 


[D} = 


‘ 4/M 
-1/M 


1/M | fO 1 1/V 0 

-7/12MJ [0 0 0 c/V 


(30) 


where p' = 1, ... ,4. 


The variables w\ and W 2 in equation (29) are given by 

w, = TiWTjtt) 

w 2 = cosA(Ti^Tjlj }Z +TiT M („ ) - 3inA(Ti in Tjij )y +TiT M („ ) 
where A is the sweep angle. 


(31) 


[Hij] in equation(29) is a matrix of order 2 xN where N = (m + l)(n + 1). A typical 
column of the matrix is given by 


f TT 1 

t cos A(Ti t , t Tj T }, x +TiTj t t£, x ) — sinA(Ti jtJ TjT) ys +T,Tj i ^, y ) J 

where k = 1, 2, . . . , N. 


(32) 


The column vector {x} in equation(29) is the vector of aerodynamic state variables 
and Pij} T is the vector of generalized displacements. 

It should be noted that in the integrations performed along the quarter-chord line in 
equation(29), a constant value of the section lift-curve slope Cn„ is used. But for a finite 
3D wing, the wing lift-curve slope depends on the planform of the wing which makes C Na 
sensitive to the shape variations of the wing. 


20 


Using equations(27),(28) and (29), equation(26) can be written as 


N 


6W nc — 'y ' QiSPi , 


t=l 


N = (m + l)(n + 1) 


( 33 ) 


where 


Qi = [C'){x} + [D 




(34) 


The aerodynamic state equations(6) for a typical section perpendicular to the quarter 
chord line were in the form 

{i} = [A]{x} + [£]{«} 

An integration of these state equations along the quarter chord line to consider the effect 
of finite span yields 




(35) 


The equation(25) can be written as a set of first order ODE’s in {P t> } and [P i} ] which 
will be represented by {p,} and {^} respectively. It can be coupled with equations(34) 
and (35) to generate the aeroelastic equations of the wing in the form 


(36) 


Since a Chebyshev polynomial of order 5 is chosen for the displacement function in tj 
and t we have 36 generalized coefficients {p,}, their 36 time derivatives { qi } and the 8 

aerodynamic state variables {x}. The stability of this system can be determined by solving 
an 80x80 eigenvalue problem. 
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5. SENSITIVITY ANALYSTS 


The aeroelastic equations obtained as a set of first order ODEs is of the form 


[P]w = [Q]w 


( 37 ) 


which could be written as 


w = [j£]w 


(38) 


where [f?] = [P] -1 [<2] 


The derivative of the ith eigenvalue with respect to the flutter speed is given by 

(4) 

dVf {e}} r {ej.} 


(39) 


where {ej} and {e' r } are the ith left and right eigenvectors respectively. is calculated 
by recomputing the [£] matrix at a slightly higher speed than V f and using the forward 
difference technique. 


Similarly, the derivative of the ith eigenvalue with respect to any parameter p 


by 


ay _ 

9p {=!F{e'} 


is given 


(40) 


can be conveniently written as 


ffl_aa± ra+ „ n -aai 

dp dp WJ + 1F1 ap 


(41) 


and can be computed analytically, where 




( 42 ) 
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The analytical derivative of the flutter speed with respect to parameter p is then given 
by 

dV f _ Real {%) 

dp Real (-§V?) (43) 

The [E] matrix is composed of mass, stiffness and aerodynamic matrices. Obtaining 
the analytical derivatives of the mass, stiffness and aerodynamic terms with respect to any 
parameter p is straightforward for the typical section. In the case of the finite wing, the 
expressions for the analytical derivatives of the [K] and [M\ matrices (equations(21) and 
(24)) are given in [35]. Since the reduced stif&iess matrix [K] is obtained from [K] by static 
condensation the analytical derivative ™ is obtained by a succession of differentiations 
using the chain rule. The derivatives of the aerodynamic terms are obtained by taking the 
analytical derivatives of those terms that are explicit functions of the shape parameters, 
given in the Appendix. 
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6. EVALUATION ANALYSIS 


The flutter characteristics of the airfoil are found by calculating the complex eigen 
values A* = <t* + iuj k at various values of free stream velocity. Flutter occurs at the lowest 
speed for which any a ^ becomes positive. 

The flutter speed was determined for the following case, the results for which have 
been presented by Leishman and Crouse [32]. The parameters used are /i = 100, x$ = 
0.25, r e = 0.5, io h = 10 rad/s n u 9 = 50 rad/s., a k = -0.5,6 = 5 m., C Na = 14.65, r ac = 
0.286, M = 0.85. Flutter was found to occur at 92.34 ft/s., i.e., a non-dimensional speed 
of V/bv 9 =4A3 which agrees well with the value of V/buo=4A reported in [32]. 

The damping ratio for each of the aeroelastic modes is given by 


(k = - 


<*k 

y/*k+ u l 


(44) 


A plot of the variation of the damping ratio ( with non-dimensional speed Vfbu> 9 is given 

in Fig.3. Fig.4 shows the variation of flutter speed predicted by this theory for different 
values of u) h /uj d . 


The flutter analysis was also carried out in the time-domain for comparison with the 
results obtained from the eigenanalysis. A time-integration of the first order ODEs repre- 
senting the aeroelastic system was done using the Wilson-0 method. The plunge and pitch 
amplitudes of motion are plotted with respect to time in Figs. 5 and 6 respectively, at three 
different non-dimensional speeds, including the flutter speed. It can be seen that at speeds 
below the flutter speed, the oscillations that axe set in due to any initial disturbance given 
to the airfoil die out as time progresses, whereas, at speeds above the flutter speed, the 
displacement amplitudes increase with time, leading to instability. At the flutter speed, 
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the oscillations are able to maintain a constant amplitude, denoting a neutrally stable 
condition. 

The sensitivity of flutter speed with various parameters namely /i, x$, r$, and u>$ 
was calculated by both analytical and finite difference methods. In the analytical method, 
the derivatives of the [E] matrix (see section on Sensitivity Analysis) with respect to the 
above mentioned parameters were calculated analytically. The finite difference derivatives 
were calculated for step sizes of 1%, 0.1% and 0.01%. The parameters were perturbed one 
at a time using these step sizes and the flutter speed recomputed. A forward diff erence 
scheme was then applied to compute the derivatives. It can be seen from the results shown 
in Table 1 that the forward difference derivatives obtained using a step size of 0.01% have 
good agreement with the analytical values. 

Figs. 7—11 show the variation of flutter speed obtained by eigenanalysis with respect 
to various parameters. In each case, the sensitivity derivative computed at a particular 
value of the parameter is also shown. 

Having gained confidence in the accurate prediction of the sensitivity derivatives, the 
sensitivity derivatives were computed both analytically and by finite difference method for 
a range of transonic Mach numbers. The C Na {M) and x ac (M) for these calculations were 
obtained from [4]. The results are tabulated in Table 2. There is good agreement between 

the analytic derivatives and the finite difference derivatives obtained using a step size of 

0 . 01 %. 

Before performing the sensitivity calculations of the flutter speed of the wing with re- 
spect to shape parameters, comparison of natural frequencies and predicted flutter speeds 
was made with results from other sources. The first three natural frequencies of an unswept 
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wing and its flutter speed in subsonic flow were compared with results reported by Lands- 
berger and Dugundji [36] for different laminate sequences in Table 3 and Table 4. The ma- 
terial used for the wing is Hercules Graphite Epoxy (ASl/3501-6) with properties: E x = 98 
x 10 9 Pa, E 2 = 7.9 x 10 9 Pa, u 12 = 0.28, G 12 = 5.6 x 10 9 Pa and p = 1520 kg/m 3 . The 
thickness of each ply is 0.134 x 10~ 3 m. The flutter data used for comparison [36] are 
the experimental results from the wind tunnel tests performed in the MIT Acoustic wind 
tunnel. The results agree fairly well. 

An experimental investigation of the flutter characteristics of the wing in transonic 
flow was performed in the Langley transonic dynamics tunnel, the results of which have 
been reported by Yates [37]. The 2.5 foot weakened 3 model is used for comparison 
with our results. The material properties of the laminated mahogany wing are: Ei = 
0.47072 x 10 6 psi, E 2 = 0.01883 x 10 6 psi, u 12 = 0.28, G 12 = 0.05975 x 10 7 psi and 
p = 0.60267 slug /ft 3 . The wing dimensions are: Area — 3.782 ft 2 , Aspect ratio = 
1.6525, Taper ratio = 0.6576, Sweep( 1/4 chord ) = 45° and thickness = 0.056 ft. The 
natural frequencies for the wing obtained from the present analysis is compared with the 
measured natural frequencies from [37] in Table 5. A good agreement is obtained. The 
flutter speed obtained from this analysis is compared with the non-dimensional flutter data 
measured in air for the WEAK3 model [37] at different Mach numbers in Table 6. The 
non-dimensionalizing parameter used is b a u> ay /p where p is the mass ratio, b s is the root 
semichord, and u) a is the angular frequency of the first torsion mode. For this model, b 3 
and u a are 0.9165 ft and 230.9 rad/s respectively. It should be noted that the 2.5 foot 
WEAK3 model has a low aspect ratio of 1.6525. In our aerodynamic model, we have 
assumed the lift forces to be acting along the quarter chord line and the moments about 
the quarter chord line, which is a good assumption for high aspect ratio wings. In reality 
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however, the aerodynamic forces and moments are distributed over the surface of the wing. 

Yet, the flutter results obtained are encouraging. 

Having achieved good prediction of flutter speeds, sensitivity analysis of the flutter 
speed of the wmg with respect to shape parameters is carried out. The wing box is shown 
in Fig.12. The wing skins are made of 0* laminated Graphite/Epoxy (T300/1V5208) with 
the following material properties: Ey = 181 x 10» Pa, E 2 = 10.3 x 10» Pa, „ 12 = 0 .28, 

G.s = 7.17 x 10 9 Pa and p = 1600 k y/m 9 . The critical airspeed of the wing is shown in 

Fig.13 as a fraction of the quarter-chord sweep angle. As seen horn the graph, divergence 

(zero frequency flutter) instability is critical upto a sweep angle of about 16» and for 

hrgher sweep angles, the flutter mode is the unstable mode. Tables 7 and 8 give the shape 

sensitivity derivatives of the divergence speeds and flutter speeds of a wing at Mach 0.9. 

The analytical shape derivatives agree weU with the finite difference derivatives obtained 
with a stepsize of 0.01%. 

The critical speeds of the wing obtained by perturbing one shape parameter at a time 

from the baseline configuration are shown in Fig.14-21. The prediction of critical speed by 

analytical sensitivity calculations is also superposed. The sensitivity derivative obtained 

forms a tangent to the critical speed curve at the value of the shape parameter at which 
it is computed. 

In order to observe the aeroelastic phenomena in real time, one of the coefficients of 
fhe displacement flmction was perturbed and the system of equations(36) was integrated 
with .aspect to time using the Wilson-* method. The wing tip displacement is plotted as a 
function of time at different speeds in Fig.22-28. The tip displacement of the unswept wing 
at the divergence sp«d is shown in Fig.22. It can be seen that the displacement approaches 
a constant amplitude at this speed. Above this speed, the displacement increases with time 
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as shown in Fig. 23. For the 15° swept wing, the displacement slowly approaches a constant 
value at the divergence speed as shown in Fig.24. The oscillatory nature is due to the fact 
that the flutter and divergence speeds are close enough. Fig.25 shows the flutter condition 
for the 15° swept wing. Since the wing has already diverged at a lower speed, we see the 
constant amplitude oscillations about a diverging mean position. Fig. 26, 27 and 28 show 
the tip displacement of the 30° swept wing below the flutter speed, at the flutter speed 
and above the flutter speed, respectively. 


The lift-curve slope of a finite wing depends on the planform of the wing and hence is 
sensitive to shape variations of the wing. For a finite wing in compressible flow, the wing 
lift-curve slope as given by Hauptman and Miloh [38] can be written as 


C Na = 


E a (h) 

k+(arc$in h)/h 


T — M 2 cos 2 A 


where k — 4/7 rAR and h = ^/(l — & 2 ) and AR is the aspect ratio of the wing. E(h ) is the 
complete elliptic integral of the second kind and is given by 



(1 — h 2 sin 2 <j>) 1 / 2 d<f> 


(46) 


Using the above expression for the wing lift-curve slope in equation(29), the flutter 
calculations were carried out. The critical airspeed of the wing is shown in Fig. 29 by varying 
the sweep angle of the wing. The sensitivity of the flutter speed of the wing with respect to 
shape parameters and the flutter speeds obtained from reanalysis by perturbing one shape 
parameter at a time are shown in Fig. 30-33. The sensitivity derivative forms a tangent to 
the flutter speed curve at the baseline configuration and gives a linear approximation to 
the flutter speeds within a certain range, about the baseline value. 
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7. CONCLUDING REMARKS 


The dynamic aeroelastic behavior of a two degree of freedom airfoil, with plunging and 
pitching degrees of freedom, was examined using a state-space approach in the transonic 
flow regime. The mdicial response functions presented by Leishman and Nguyen [14] were 
used to represent the unsteady aerodynamic behavior. The flutter speed of the airfoil 
was determined using an eigen analysis of a system of 12 first order ordinary differential 
equations. The results obtained from the eigenanalysis was compared wth the results 
obtained by performing a time-integration of the aeroelastic equations. The sensitivities of 
the flutter speed with respect to the mass and stiffness parameters were computed by both 

analytical and finite difference methods and the results obtained are in excellent agreement 
with each other. 

The compressible unsteady aerodynamic theory using indicial response functions was 
successfully modified to represent the aerodynamic forces and moments on a finite span 
wing. Using this aerodynamic state-space model and the structural formulation based on 
Ritz technique, flutter analysis of wings were carried out in transonic flow. The use of 
Chebyshev polynomials for Ritz functions gives the added benefit of closed form analytical 
expressions for the derivatives of stiffness and mass matrices with respect to the shape de- 
sign parameters of the wing. This avoids the uncertainities and the computational expense 
associated with finite difference derivative calculations. The shape sensitivity derivatives 
of the critical speed of the wing were computed by analytical and finite difference meth- 
ods and they are in excellent agreement with each other. These shape derivatives of the 
flutter response of a wing in transonic flow would be very useful to a designer in the initial 

design phase, thus avoiding the necessity of a reanalysis for small changes in the design 
parameters. 
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APPENDIX 


The aij's are given by 


where 


2V 

a n = 

c 

, 2 F 2 , 

a 22 = -( — )jS 2 b 2 

c 


a 33 = 


a 44 = 


K a Tj 

1 


K q Tj 

a 55 = ~ (^3^aAfT/) -1 

066 = -{WK a MTj)~ l 

077 = -bs/3 2 ( — ) 

c 


088 = ~ 


Kq M Tl 


-1 


Ka{M) = [(1 - M) + TT/3M 2 (A 1 b 1 + A 2 & 2 )] 
K q (M) = [(1 - M) + 2ir0M 2 (Aibi + A^)]- 

Kq\((M) = 


15(1 -M) + 3irflM 2 bs 


The c./s axe given by 


2V 

c u = C Na (—)0 2 A 1 b 1 

2V 

c i2 = C Na (-—-)0 2 A 2 b 2 


4 , 

c “ = 77(- 


-) 


K a Ti 

= M^~K q T^ 

2V 

c 2i = C Na (~) i 3 2 A 1 b 1 ( 0.25 - x ac ) 
2V 

c 22 = C Na (~— )£ 2 4 2 6 2 (0.25 - x ac ) 
•^3055 

M 

■^ 4^66 


C 2 5 = ~ 
c 26 = 
c 27 = 


M 

CN a 


u ^ 2 (f> 

C 28 — TnTT? (“ ~ ) 


UM y K qM T! 

The constants are given by = 0.3, 4 2 = 0.7, 4 3 = 1.5,^ 
0.53, 63 = 0.25, 6 4 = 0.1, 65 = 0.5. 


(49) 


= —0.5, 61 = 0.14,62 = 
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Analytical derivatives 


Aspect ratio (AR) : Area (5), Taper ratio (<r), Sweep (A) 


The wing coordinates are aq, x 2 , £3 and X 4 . 



x\ — 0.75 cr 

x 2 = span tan A + 0.75 ct 
X 3 = span tanA — 0.25 ct 
X 4 = — 0.25 cr 
p = (#1 + x 2 + #3 + x a) 
PP = (x 2 + X 3 ) - (x\ + £ 4 ) 
crct = cr + ct 
rt = cr — ct 


For any point on the quarter chord line ( — 1 < rj> < 1 ) 

y = 0.5 span (1 + xl>) 
x = y tanA 



span 

(£ pp — 4z + p) 

77 = — - 

(£ rt — crct ) 

c = [cr + (ct — cr)(l + V’) 0.5] cosyl 
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Local chord (c) 


_dc_ = _ S°- 5 [l - (1 - tr)(i + 0) 0.5] cos A 
dAR 1^5(1 + tr) 

— = SlzSLz £Xj + 0) 0-5] co^yl 
dS 7M 5(1 + tr) 

dc cr 0 cos/1 
5<r (1 + tr) 

dc 

-QA = — c £an/l 
&= ( g i +X-7+Z, +xA 


dp 

_ 0.5 5 0 5 r 

dAR 

AR°* 

dp 

0.5 p 

dS 

S 

dp 

no 

dtr 


dp 

2 span 

dA 

cos 2 A 


tanA — 


AR 


gP — fx-7 + g’t'l — fXl + X,«) 


dpp 

dAR 

d pp 

dS 


= 0.5 S 05 
0.5 pp 


2 tanA 


+ 


(1 - tr) 


AR ' 0 5 (1 + tr) AR 1 


.5 


dpp — cr 
dtr (1 + tr) 
dpp _ 2 span 
dA cos 2 A 


crct = cr 4- r.t 


dcrct S 05 


dAR AR 1 - 5 

dcrct 0.5 crct 


dS 

dcrct 

dtr 

dcrct 

~dT 


= 0.0 
= 0.0 
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rt = cr — ct 


drt _ 5° 5 (1 - fr) 

“ AR 1 -* (1 + fr) 

drt _ 0.5 rt 
55 ~ ~S~ 


drt _ 2 cr 

dtr (1+tr) 



x_=JL tanA 


dx _ S 0,5 tan./l(l 4- 0) 

&4.R 4 >1R 0 - 5 

dx 

-QS = 0.5 (1 + %!)) span tanA 
dx 

^ = 0.0 

atr 

dx _ 0.5 span (1 + 
d/1 cos 2 /l 


?? = (£ pp - 4s + p}/(£ rt - crct) 


For any parameter v, 


drj 

dv 

dZ 

dv 


((_rt crct)(£^g -4 jf + |g) - (£ pp - 4a + p)(f jg 


dcrct ^ 
JW 
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Table 1. Sensitivity of flutter speed with respect to various parameters (M=0.8S) 

(/i — 100, x# — 0.25, r$ = 0.5, = 10 rad/s., 

ojq = 50 rad/s. , 6 = 5 in., a h = -0.5) 


Parameter 

Analytic 

derivative 

Finite Difference 
derivative 

1%“ 

0.1% 

0.01% 


0.39296 

0.38798 

0.39244 

0.39289 


-253.60673 

-248.75989 

-253.10277 

-253.54363 

ro 

239.41460 

236.53840 

239.11255 

239.37244 

Uh 

-4.35016 

-4.31603 

-4.34651 

-4.34958 

a ; — l: j. 

2.71689 

2.68798 

2.71383 

2.71644 


a indicates step size 




Table 2. Sensitivity of flutter speed with respect to various parameters 
at different transonic Mach numbers 
(/i = 100, x$ = 0.25, r$ = 0.5, = 10 rad/s., 

u>e = 50 rad/s., b = 5 in., a h = -0.5) 


Parameter 

Mach 

number 

Sensitivity derivatives 

Analytic 

Finite Difference 


0.1%* 

0.01% 


0.8 

0.44303 

0.44247 

0.44298 


0.85 

0.39296 

0.39244 

0.39289 


0.8625 

0.46202 

0.46143 

0.46196 


0.87 

0.53132 

0.53010 

0.53070 

xg 

0.8 

-241.52263 

-241.07122 

-241.47758 


0.85 

-253.60673 

-253.10277 

-253.54363 


0.8625 

-463.70394 

-462.54019 

-463.58734 


0.87 

-734.98634 

-731.92101 

-733.99406 

r$ 

0.8 

243.06012 

242.76615 

243.03092 


0.85 

239.41460 

239.11255 

239.37244 


0.8625 

308.46665 

308.08430 

308.42838 


0.87 

391.63750 

390.75475 

391.18409 


0.8 

-3.78605 

-3.78311 

-3.78576 


0.85 

-4.35016 

-4.34651 

-4.34958 


0.8625 

-6.90676 

-6.90064 

-6.90615 


0.87 

-10.17027 

-10.14904 

-10.15867 


0.8 

2.80282 

2.79985 

2.80253 


0.85 

2.71689 

2.71383 

2.71644 


0.8625 

3.58044 

3.57671 

3.58007 

b j . 

0.87 

4.60379 

4.59453 

4.59857 


b denotes step size 



Table 3. Comparison of natural frequencies of an unswept wing 
(Area = 0.02318 m 2 , Aspect ratio = 4.0132, Taper ratio = 1.0) 
for different laminate sequences 


Laminate 

sequence 

Natural frequencies (Hz.) 

Present 

Landsberger and Dugundji [36] 

First 

Second 

Third 

First 

Second 

Third 

[0 2 /90], 

11.03 

39.30 

69.06 

10.8 

39 

67 

[15 2 /0], 

8.86 

42.62 

63.25 

8.5 

48 

58 

[±15/0], 

10.12 

48.9 

64.94 

9.9 

50 

63 

[+30 2 /0], 

6.21 

37.57 

57.78 

6.0 

41 

60 

[±30/0], 

7.73 

48.76 

64.42 

7.8 

50 

65 


Table 4. Comparison of flutter speed of an unswept wing 

(Area = 0.02318 m 2 , Aspect ratio = 4.0132, Taper ratio = 1.0) 

in subsonic flow 


Laminate 

sequence 

Flutter speed (m/s) 

Present 

Experimented [36] 

[0 2 /90], 

24.9 

26 

[15 2 /0], 

23.2 

25 

[±15/0], 

28.1 

28 

[±15/0], 

18.9 

21 

[+302/0], 

28.8 

29 




Table 5. Comparison of natural frequencies of the 2.5 foot WEAK3 model 

(Area = 3.782 ft 2 , Aspect ratio = 1.6525, Taper ratio = 0.6576, Sweep = 45°) 

with measured values 


Natural 

frequencies (Hz.) 

Present 

Experimental [37] 

First 

9.91 

9.60 

Second 

36.75 

38.10 

Third 

60.78 

50.70 

Fourth 

102.11 

98.50 


Table 6. Comparison of flutter speed of the 2.5 foot WEAK3 model 
(Area = 3.782 ft 2 , Aspect ratio = 1.6525, Taper ratio = 0.6576, Sweep = 45°) 
with flutter data measured in the transonic wind tunnel 


Mach 

Mass 

Density 

Flutter 

Non-dimensional 

number 

M 

ratio 

of air 
(slugs/ ft 3 ) 

speed, V 
(ft/s) 

speed (V/b 3 w Qy /fT) 



Present 

Experimental [37] 

0.499 

f 33.465 

F 0.000830 

476.90 

0.3896 

0.4459 

0.678 

68.753 

0.000404 

640.00 

0.3647 

0.4174 

0.901 

143.92 

0.000193 

829.44 

0.3267 

0.3700 




Table 7. Sensitivity of divergence speed of the wing at M=0.9 
(Area = 20 m 2 , Aspect ratio = 10, Taper ratio = 0.5) 
with respect to shape parameters 


Sweep 
angle (deg) 

Divergence 

speed(m/s) 

Parameter 

* * 

Analytic 

derivative 

j™ Finite difference derivative 

1.0% a 

0.1% 

0.01% 

0 

137.41 

Aspect ratio 

-4.9672 

-4.9276 

-4.9630 

-4.9642 



Area 

-4.2460 

-4.2156 

-4.2418 

-4.2433 



Taper ratio 

7.2069 

7.4086 

7.2601 

7.2308 



Sweep angle 

147.0488 

147.1087 

147.2050 

148.5662 

15 

203.74 

Aspect ratio 

16.9889 

17.4989 

17.0278 

16.9773 



Area 

-18.3916 

-17.9050 

-18.3307 

-18.3768 



Taper ratio 

440.3419 

451.3756 

441.1370 

440.0515 



Sweep angle 

434.5355 

437.5666 

434.5507 

434.0796 


a indicates step size 




Table 8. Sensitivity of flutter speed of the wing at M=0.9 
(Area = 20 m 2 , Aspect ratio = 10, Taper ratio = 0.5) 
with respect to shape parameters 


Sweep 


Flutter 

Parameter 

Analytic 

Finite difference derivat: 

speed(m/s) 


derivative 

0.1%“ 

0.01% 

229.90 

Aspect ratio 

-12.6812 

-12.6675 

-12.6779 


Area 

-6.4607 

-6.4563 

-6.4592 


Taper ratio 

-184.8167 

-184.7288 

-184.7691 


Sweep angle 

-104.4327 

-104.3587 

-104.4143 

209.67 

Aspect ratio 

-10.1067 

-10.0929 

-10.1005 


Area 

-4.3196 

-4.3181 

-4.3170 


Taper ratio 

-199.1987 

-199.1011 

-199.0925 


Sweep angle 

-37.7435 

-37.6769 

-37.5526 

213.22 

Aspect ratio 

-6.6440 

-6.6311 

-6.6404 


Area 

-5.0568 

-5.0542 

-5.0554 


Taper ratio 

-172.6965 

-172.6105 

-172.6414 


Sweep angle 

67.8415 

67.9632 

67.8975 


1 indicates step size 




Fig 1. Two degree of freedom airfoil 
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Fig. 5a Plimge displacement Vs time 
(1 time step = 0.01 sec.) 
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Fig. 5b Plunge displacement Vs time 
(1 time step = 0.01 sec.) 
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Fig. 5c Plunge displacement Vs time 
(1 time step = 0.01 sec.) 




(1 time step = 0.01 sec) 
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Fig. 6c Pitch displacement Vs time 
(1 time step = 0.01 sec.) 
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Fig. 9 Flutter speed Vs Radius of gyration 
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Fig 12. Wing box section used 
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g. 13 Critical speed Vs Sweep Ang 

(AR=10, Area=20 m, TR=0.5) 
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Fig 14. Divergence speed Vs Aspect ratio (M=0.9) 

(AR=10, Area=20rrr, TR=0.5, Sweep=OP) 
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Fig 15. Divergence speed Vs Area (M=0.9) 

(AR=10, Area=20m 2 , TR=0.5, Sweep=0°) 
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Fig 16. Divergence speed Vs Taper ratio (M=0.9) 

(AR=10, Area=20m 2 , TR=0.5, Sweep=CP) 
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Fig 17. Divergence speed Vs Sweep angle (M=0.9) 

(AR=10, Area=20m 2 , 777=0.5, Sweep=(f) 
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Fig 18. Flutter speed Vs Aspect ratio (M=0.9) 

(AR=10, Area=20rrf, TR=0.5, Sweep=15°) 
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Fig 21. Flutter speed Vs Sweep angle (M=0.9) 

(AR=10, Area=20nr, TR=0.5, Sweep=1&>) 
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Fig 24. Displacement Vs Time at V=203.74 m/s 

(AR=10, Area=20m 2 , TR=0.5, Sweeps 15 °) 
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Fig 25. Displacement Vs Time at V=209.67 m/s 

( AR=10,Area=20m ?, TR-0.5, Sweeps 15°) 
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Fig 27. Displacement Vs Time at V=213.22 m/s 

(AR=10, Area=20m 1 , TR-0.5, Sweep-30 0 ) 
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Fig 28. Displacement Vs Time at V=220 m/s 

(AR=10, Area=20m 2 , TR=0.5, Sweep=30°) 
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Fig. 29 Critical speed Vs Sweep Angle 

(AR=10, Area=20 m, TR=0.5) 
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Fig 31. Flutter speed Vs Area (M=0.9) 

(AR=10, Area=20m 2 , TR=0.5, Sweep=1&) 
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Fig 32. Flutter speed Vs Taper ratio (M=0.9) 

(AR=10, Area=20m 2 , TR=0.5, Sweep=lF) 



Flutter speed (m/s) 



Sweep angle (degrees) 

Fig 33. Flutter speed Vs Sweep angle (M=0.9) 

(AR=10, Area=20m 2 , TR=0.5, Sweep=15°) 



